Decoherence and purity of a driven solid-state 
qubit in Ohmic bath 

Zhi-Qiang Li, Xian-Ting Liang* 
Department of Physics and Institute of Modern Physics, 
Ningbo University, Ningbo 315211, China 

March 31, 2008 



Abstract 

In this paper we study the decoherence and purity of a driven solid- 
state qubit in the Ohmic bath by using the method based on the master 
equation. At first, instead of solving the master equation we investigate 
the coefficients of the equation which describe the shift in frequency, dif- 
fusive, decoherence, and so on. It is shown that one of the coefficients (we 
called it decoherence coefficient) is crucial to the decoherence of the qubit 
in the model. Then we investigate the evolution of the purity of the state 
in the model. From the analysis of the purity we see that the decoherence 
time of the qubit decrease with the increase of the amplitude of the driven 
fields and it is increase with the increase of the frequency of the driven 
fields. 
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1 Introduction 

Solid-state qubits are considered to be promising candidates for realizing build- 
ing blocks of quantum information processors. In past few years, many kinds 
of these qubit models are introduced and many efforts not only theoretical but 
also experimental have been performed on investigating the decoherence, relax- 
ation and manipulation for these qubit models. Most of the investigations for 
these models are considered them as a spin-boson model which is modeled with 
a two-level system coupled to a bath and the bath always constructed with a 
set of hamonic oscillators. The qubits can be manipulated with a driven field 
[U [2 [3l |H [5] , so a more realistic description requires the inclusion of an external 
control field in the qubit-bath model. The quantum tunneling or other quantum 
properties of the driven two-level system have been investigated at dept in last 
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years [5J [7] . Recently, the driven spin-boson model attracting a lot of attention 
because of its interest in connection with quantum computing with solid state 
devices. The control of the coherence dynamics of a two-level atom placed in a 
lossy cavity with an external periodic driving field [21] , the dynamics of a XOR 
gate operation with an external source [25], and the consequence of driving in 
terms of multi-phoyo transtions experimentally |26| are investigated. In this pa- 
per, we shall analyze the influence of the driven control field on the decoherence 
of the qubits by using the master equation method [I2J Q21 [THl HH] ■ At first we 
shall derive out a master equation from the Hamiltonian of the driven qubit in 
the Ohmic bath. Then we analyze the decoherence and other decay coefficients 
based on the master equations. Finally, a rate of purity decay for this model is 
obtained and discussed. 



2 Models and Master Equation 

In this paper we concentrate on the case of a persistent current qubit based on 
Josephson junction [9l [10l [11] . It is driven with the magnetic flux through the 
loop and damped predominantly by flux noise with Gaussian statistics. This 
setup is accurately expressed by the driven spin-boson (DSB) model [14] , fT6| . 

m 

H = H S + H B + Hj, (1) 



where 



H s = \e{t)a z - ±Aa Xl 
h b = + m n u 2 n x 2 n ),Hi = a z J2 n A„av, 



(2) 



Here, Oi (i = x, y, z) are Pauli matrices. The quantum environment is modeled 
with an infinite set of harmonic oscillators of mass m n , angular frequency u) n , 
momentum p n and position coordinate x n which are coupled independently to 
the spin a z with strength measured by the set {A„} and e{t) — e + scos(fW) 
is the external, time-dependent control field with the static basis sq. For s = 0, 
the system degenerates to be a common two-level one. 

The evolution of the total density matrix for the system, in the interaction 
picture, reads 

ihp T = [Hi,p T ], (3) 
where the interaction representation of the operators are given by 

~ /.x ,iH t -iH Q t .. 

p T (t) = exp(— — )p T exp( — - — ),o 

Hl (t) = exp(^)i? / exp(^-^). (4) 

Here p T = p T (0) and Hi = Hi(0), where Hi (0) is the Hamiltonian of the 
system in the Schrodinger picture. The perturbative expression of the Eq. |4]) 
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in the interaction picture is given to second order of Hi by 



p T (t) = p ~(o) + i/ dh&ih),^)} 

+ {\? I dh I' ' dtzlHjfrUH^firiO)]]. (5) 
w Jo Jo 

The reduced density operator for the system is defined by p(t) ~Trg(p T (t)) : 
Where Tr^ indicates a trace over environment variables, then we assume that 
initially the system and environment are uncorrelated that the total density 
matrix is a tensor product of the form p T (t) = p(0) <£> Pb(0). So we get 
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in Jo 

+ (h 2 I dt x r dt^Hjih), [Hj(t 2 ), p(0) <g> p B (0)]]. (6) 

''' Jr. Jo 

Next, we make a rather trivial operation that enables us to finish the derivation 
in a simple way, we can express the initial state p(0) in terms of p{t) using the 
same perturbative expansion and rewrite the Eq. ([5]) while the initial state p(0) 
appears in the right-hand side only, then inserting this expression into Eq. |6j) 
and making the derivation [20] . we can obtain 

p(t) = 4/ &&[Ei(h),P ® Pal ~ I dt 1 Tr[Hi(t),[H I (t 1 ) ) p®p~ B }} 
*«■ Jo n Jo 

+ _L [ dhTr&it), TM#i(*i), P ® Pb\ ® PbI (7) 

n Jo 

The Eq. ([7]) is in the interaction picture. By virtue of the evolution operator 
U s — exp[— j- dt\H s {ti)\ we rewrite the equation in the Schrodinger picture, 
and considering the coupling term in the DSB with Hi — a z J2 n ^nX n , then we 
get the master equation of DSB in the Schrodinger picture is 



= —[H e f f ,p] + ^X^[ A " ( x n) a z,P\ 
n 

1 /"* 

"12 E / dt 1 (k^[a z ,[a z (t 1 -t),p]]+k^[a z ,{a z (t 1 -t),p}]),(8) 
n „ Jo 



with 



k { n ] = ^A* ({x n (t), x n (ti)}) - X 2 n (x n ) (x n ) , 

fci 2) = \>?n ([Xn(t),X n (h)}) 

H e j j = a- x a x + b x Oy + C\Sr z , 



fci 2) = \x 2 n ([x n (t),x n (h)]), (9) 
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where, 



A e t + ^sin(fti) 1 



01 = -cos( ^ ) + — A(s t + ^ sin(fti)) 



s 



ft 



+Ai(e + s cos(fti))] x cos( — ) sin( 
A . .s t + ^sin(fti) 1 



At N _. ,e t + ^sin(fti), 



ft y v ft 



61 = — sin( ^ — A(e i + -sin(ftt)) 



4ft 



s 



ft 



+At(e + a cos(fti))] x cos(^) CQg( ^ + yn(ftt) ^ 



A, 1 



ci = - — (e + scos(ftt)) + — A(e t + — sin(fti)) 



2" v " 4ft 



s 



ft 



+At(e + scos(ftt))]sin(^). (10) 

Here, the notation [•,•], {•,•} denote the commutators and anticommutators, 
and (..) =Tr(..p) are reserved for quantum expectation values. So far, we only 
make two important assumptions: first, we used a perturbative expansion up 
to second order, which is accurate enough as the interaction small enough com- 
paring to the Hamiltonians of system and bath; second, we assumed that the 
initial state is not entangled. Next we consider another important approxima- 
tion, which is usually considered in the master equation, the Markov approxi- 
mation. Assuming the kernel kn are strongly peaked at the point t = ti, with 
slowly varying functions, then we can transform the temporal integrals over the 
variable r = t — t\. And assumed that the initial state of the environment to 
be thermal equilibrium at temperature T = l/ksP, thus the first order term in 
the master equation of DSB disappears because (x n ) — 0. Therefore, the master 
equation becomes 



1 r„ , 1 



P = -^[Heff,p] --gj dti{y(t{)[<j z ,[<j z {-ti),p]] - w7(ii)[<7z,{<T 2 (-ti),p}]), 

(11) 

with the two kernels are the dissipation and noise kernels respectively and are 
defined as 

= ^Z! A n({*n(*).Zn(0)}> = jf du J(w) COs(ijt) COth(^), (12) 
1 f 00 

1® = 2h S A ™ tt Xn ^> Xn( - ^ = / dujJ ^ sin ^' (13) 

n ^ 

where J(oj)is the bath spectral density function defined by 

J H = ?E^^--«)- ( 14 ) 
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We consider that the Hamiltonian of the spin system H s = ^e(t)cr z — \A& X , 
the evolution operator has the form 



U s = cxp[-- J dhHsih)] 



IS % 

= exp{- — [e Q t+ -s'm(nt))a z }exp{— Ata x } 

X 6Xp{ " Wf [£ot2 + *S sin{nt ^v + °(* 3 )}' (15) 

where we using the formula e A+B = e A e B exp^V E 2 jlV.jmri[ B ^ Aii) ^ BU) }- 
If we set the time t very small (namely, t is smaller than the characteristic time 
of the qubit), we can ignore the term of o(t 3 ). So we can solve the Heisenberg 
equations for the system and determine the operator a z (t) to be 



.At. .At , s fn N1 , ,Ai, 
a z {t) = -a x cos{— )sin(^[e t + — sin(Ot)J) - er y sm( — ) 

+<j z cos(-^) cos(-^[eoi + 75 sin(fit)]). (16) 



So the results of this paper are all the short-time results. Substituting this 
equation into Eq. (|lip . we obtain the final expression for the master equation 
of DSB. It is 

P = ^[H eff + n(t),p]-D(t)[a z ,[a z ,p]]-G(t)[a z ,[a x ,p]} 

-fitflvz, Wy,p]] +iri(t)[<T z ,{<J x ,p}}+ir 2 (t)[<T z ,{<Ty,p}} (17) 



where 



D(t) 



J dti^(ti)(cos(^i)cos(^i[e ti + ^ sin(fiti)])), 



./(/!= / dt l! /(t 1 )(- si n(^)), 



(V!/! / )(-cos(^)sin(^[£ t + ^sin(fi/)])). 



fi(t) = / dtir7(ti)(cos(^i)cos(^i[e ti + ^sin(Ot 1 )])), (18) 



n(t) = / dti77(ti)(-cos(-^)sin(^[e t + ^sin(Ot)])), 



r 2 (t) = / rfti»?(t 1 )(-sm(— )). 
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3 Decoherence and purity decay 



Having obtained a master equation describing the evolution of a controllable 
solid state qubit system coupled to an Ohmic bath, we now proceed to examine 
the consequences of that evolution. In the following we shall analyze the problem 
with two different methods. At first, we use the coefficients in the master 
equation to detect the instantaneous effects of the environment. Then we can 
investigate the time evolution of the purity. For the first method, all the effects 
including controlled external field and uncontrolled environment are considered 
in the above coefficients. These coefficients are to renormalize the frequency 
as well as to introduce the decay of the system. From these equations, it is 
possible to have a qualitative idea of the effects for the environment producing 
on the system. First, we can observe that the effective Hamiltonian H e ff is 
evolved with time-dependence. The system Hamiltonian Hs and the coupling 
Hamiltonian Hi do not commute with each other and the evolution of H e f f is 
periodic because the driven controlled field is periodic. The term fi(i) is the 
shift in frequency which produces the renomalized frequency. This term does not 
affect the unitarity of the evolution. The terms D(t), f(t), G(t), r\(t) and r2(t) 
are diffusive terms and bring about non-unitary effects. The terms of r\{t) and 
r2(i) are the dissipation coefficients related to the dissipation kernel u(t) defined 
already, which play the role of a time-dependent relaxation rate, are independent 
of temperature. And D(t), f(t) and G{t) are the diffusion coefficients, which 
produce the decoherence effects, they are not only time-dependent but depend 
on temperature. Of cause, the explicit time dependence of the coefficients can 
only be computed once we specify the spectral density of the environment. To 
illustrate their qualitative behavior, we focus on the case that the environment 
is the Ohmic bath in the following, whose spectral density of the Ohmic bath 
can be expressed as 

J(uj) = 27rawexp(-^), (19) 

where A is the physical high-frequency cutoff, which represents the highest fre- 
quency presented in the environment and the parameter a is dimensionless 
parameter reflecting the strength of dissipation. We set a = 0.01 as in [23]. We 
can plot the diffusive terms D(t), f(t) and G(t) as Fig.l. 

Fig.l 

It is shown that all these plots appear to be a periodically diverging function, 
which stems not only from the periodic control driven field but also from the 
bath, and the peaks of D(t) is much higher than that of f(t) and G(t). It 
means that the term proportional to D(t) in Eq. (fT5|) plays the main role to 
decoherence. The form of the decoherence coefficient D(i) [2D] is not necessarily 
intuitive, so we evaluate it numerically as Fig. 2. 

Fig.l 

One of the parameters (s, T, fi) is varied in Figs. 2 (a, b, c). All of the plots 
demonstrate the same basic behavior: D(t) decreases with periodic oscillations. 
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It is shown that the influence of the heat bath to the decoherence of the 
qubit can be described by these coefficients of the master equation. It can also 
be described by a computationally convenient way, namely the linear entropy 
C = tr(p — p 2 ) which is another measure of the purity of a quantum state. 
For the pure state it is approximative equals to C = 1 — MP 2 ) so we shall 
study the evolution of the purity of the system as measured by £ = tr(p 2 ) for 
simplicity [201 HZl HI] • It is equal to one for a pure state and decreases when the 
state of the system gets mixed because the destruction of quantum coherence is 
generated by evolution. By virtue of the master equation, we can easily obtain 
an evolution for the purity £ and the equation we obtained is: 

£ = -AD(t)Tr{p 2 a 2 z - pa z pa z ) - AG{t)Tr{pa x pa z ) - f(t)Tr(pa y pa z ) 

-Ar 1 {t)Tr{p 2 a y ) + A ri {t)Tr(p 2 a x ). (20) 

Setting the initial state is pure then p 2 — p. Similar to Ref. [55] we can obtain 

£(t)=~D(t). (21) 

It is interesting that the Eq. ([21]) only correspond to the decoherence coefficient 
D(t). It is shown that the purity decay rate is proportional to the decoherence 
coefficient D(t), and do not include any other coefficients. It confirms that D(t) 
plays the main role in decoherence which is analyzed in the above section in 
detail. From Eq. (f2Tj) we can also see that the larger the rate of the purity is, 
the smaller the decoherence becomes. In Fig. 3, we plot the evolutions of the 
state purity with time t. 

fig-3 

It is show that the purity of qubit decreases fast with the increasing of the 
amplitude and decreasing of the frequency of the driven control field. This is 
agree with the result in Ref. [55J that the low frequency driven field is destructive 
to the coherence of the qubit. 

4 Conclusions 

In this paper we have investigated the effects of the environment for a solid state 
qubit coupled with a driven control field. At first, instead of solving the master 
equation we investigated the coefficients in the equation which express the shift 
in frequency, diffusive, decoherence and so on. Here we suppose the environment 
is modeled with a Ohmic bath and it is coupled with the qubit linearly. It shows 
that the term proportional to the coefficient D{t) in master equation plays the 
main role to the decoherence. Then, we concretely investigated the decay of 
purity. It is shown that the decay rate of purity behaves fast as the amplitude 
increasing and the frequency decreasing of the driven field. It is suggested that 
if the driven field is necessary for some qubit system a higher frequency may 
help to decrease the decoherence. 
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5 Figures Captions 

Fig. 1: Coefficients in master equation versus time t. (a) D(t); (b) f(t); (c) 
h{t). Here, the parameters are eq — A, s — A, fl — A, A = 36A, A = 5, T = 30 
mK, here and in the following figures the times are expressed in unit of 5/A. 

Fig. 2: Dependence of D(t) on coefficients of the driven field and temperature 
T. The parameters s, 0, (3, and A are same as Fig. 1. (a) D{t) versus s; (b) 
D(t) versus 0; (c) D(t) versus = 1/kT. 

Fig. 3: Dependence of purity of the system on coefficients of the driven field. 
The parameters s, ft, (3, and A are same as Fig. 1. (a) £ versus s; (b) £ versus 
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